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ABSTRACT 

We present the mode identification of frequencies found in spectroscopic observa- 
tions of the 7 Doradus star HD 135825. Four frequencies were successfuhy iden- 
tified: 1.3150 ± 0.0003 d-i; 0.2902 ± 0.0004 d"!; 1.4045 ± 0.0005 d"!; and 
1.8829 ± 0.0005 d'^ These correspond to {I, m) modes of (1,1), (2,-2), (4,0) and 
(1,1) respectively. Additional frequencies were found but they were below the signal- 
to-noise limit of the Fourier spectrum and not suitable for mode identification. The 
rotational axis inclination and vsvai of the star were determined to be 87° (nearly 
edge-on) and 39.7 kms~^ (moderate for 7 Doradus stars) respectively. A simultaneous 
fit of these four modes to the line profile variations in the data gives a reduced x^ of 
12.7. We confirm, based on the frequencies found, that HD 135825 is a bona fide 7 
Doradus star. 



1 INTRODUCTION 



The pulsation of a star is reflective of its internal structure. 
Since the interior of a star cannot be directly seen, the anal- 
ysis of the pulsation is the most advanced way to determine 
interior properties and dynamics. Because pulsations occur 
in many different types of stars, in many stages of stellar 
lifetimes, asteroseismology can be used to further our un- 
derstanding of stellar evolution. 

Asteroseismology is the study of vibrational physics in 
stars, specifically pulsational behaviour in order to derive 
interior structure properties. The characteristics of pulsation 
are described by the "quantum numbers" n, I and m. The 
radial order, n, describes the number of interior nodal shells, 
I gives the number of nodal lines on the stellar surface and 
the azimuthal number, m, defines the number of nodal lines 
that intersect the pole. 

By analysing the non-radial pulsations and determin- 
ing the spherical harmonic pulsation modes, we obtain as- 
teroseismic information about the deep interior regions of 
stars. Non-radial pulsations are classified by the restoring 
force of the motion: gravity (g-mode) or pressure (p-mode). 
Gravity-mode pulsations propagate as deep as the convec- 
tive core interface making them ideal probes of core physics. 
The detailed frequency spectrum of excited modes places 
severe constraints on the physical conditions within the re- 
gions where these modes propagate. Successful mode iden- 
tifications can be found in [Aerts et al.| ( |2004[ ), |Zima et al.| 
1 2006 1 and Bohm et al. ( 2009 1 for various types of pulsating 



star. 



ing feature of the 7 Dor class is the presence of high-order 
non-radial g-mode pulsations in an A-F type star that is on 



or near the main sequence (see Kaye et al.||1999 and Kaye 
2007 Pollard|2009 for reviews). Thus these stars are slightly 



hotter and slightly larger than the Sun. A typical feature 
of the 7 Dor class is the presence of multiple g-mode pulsa- 
tions at frequencies lower than the fundamental radial mode 
F] Observed g-mode frequencies generally fall between 0.3 - 
3.3 d"'^, although the observed frequencies are dependent 
on the stellar rotation and direction of the travelling wave. 
The pulsation frequency range for 7 Dor stars is similar to 
the g-mode pulsators of the Slowly Pulsating B (SPB) class, 
so an estimate of the Teff or spectral class is required in 
order to distinguish between these classes. Presently there 
are less than 100 bright bona fide 7 Dor stars (see list in 



Henry et al. 20111 with a further 100 candidates thus far 



2010 Uytterhoeven 



reported by Kepler (Grigahcene et al. 
et al.|2011| ). 

The 7 Dor pulsators reside in that part of the 
Hertzsprung-Russell diagram between the classical instabil- 
ity strip and the solar instability region. Thus they bor- 
der the regions of the 5 Scuti (p-mode) pulsators and the 
solar-like oscillators. A few hybrid 7 Dor/5 Scuti stars have 
now been studied (e.g. Henry & Fekel 2005 Uytterhoeven 
et al. 2008) and further candidates proposed and still be- 
ing discovered by Kepler ( [Uytterhoeven et al.|2011 Balona 
fc Dziembowski|2011 1 . Early indications show that 7 Dor/ J 



Scuti stars may challenge our understanding of the instabil- 
ity strip (Uytterhoeven et al.||2011l. 



The study presented here focuses on one g-mode pul- 
sator, the 7 Doradus candidate star HD 135825. The defin- 



Typical radial pulsation frequencies for A-F stars are 8-24 d 



Solar-like variability is driven in the convective envelope 
of a star from stochastic excitation and has short frequen- 
cies (e.g. the F-type star Procyon has solar-like variability 
freq uencies in the range 300-1400 iiHz, equivalent to 26-121 

Bedding et al.|2010 1 . Solar-like oscillations also exist in 



stars slightly hotter than the Sun ( Michel et al.|[2008 l and 
have been discovered in the S Scuti star HD 187547 ( An 



toci et al.|2011 1. This indicates that the convective envelope 
operates efficiently in stars up to 2M0. It is possible that 
in 7 Dor stars, being located between solar-like stars and 
S Scuti stars in size and temperature, solar-like oscillations 
may also be present. No solar-like/7 Dor hybrids have yet 
been discovered but it may be possibile for high-resolution 
satellite photometry to be able to detect the high frequency 
pulsations. 

The multitude of high-order, low-degree g-modes 
present in 7 Dor stars give us information about their in- 
terior structure and physical conditions. These stars have 
both a convective envelope and a convective core, in between 
which a radiative zone exists. Driving occurs by the "flux- 
blocking mechanism" where convection at the base of the 
thin convective envelope modulates flux throughput from 



the core to the stellar atmosphere ( Guzik et al. 2000 Dupret 



et al. 20041. Asteroseismic study of interior structure re- 



quires frequencies to be identified and corresponding modes 
to be fully characterised with {l,m) values. Mode analysis 
can also be used to constrain stellar parameters such as in- 
clination ( [Wright et al.|2011 l. 

High precision photometry has been successful in iden- 
tifying up to hundreds of frequencies in 7 Dor stars (see 
Chapellier et al. (2011 1 for an example of 840 identified fre- 



quencies using the COROT photometric satellite). Photom- 
etry is invaluable in the study of 7 Dor pulsations as it 
can be used to find such large numbers of frequencies with 
high precision and also determine the corresponding I values. 
Spectroscopy has thus far been limited to ground-based ob- 
servations and the signal-to-noise requirements for study of 
spectra have meant longer integration times. The benefit of 
spectroscopic observations is the determination of both I and 
m values for each pulsational frequency. As such photometric 
and spectroscopic studies of these stars are complementary. 
To date only a handful of 7 Dor stars have spectroscopic 
mode identifications published with only a few (less than 
six) frequencies each. The lack of detailed mode identifica- 
tions has previously limited interior structure studies but it 
is hoped that identifications such as in this paper can be 
incorporated and used to refine current stellar models. 

A further aspect of interest in 7 Dor stars is the effect 
of rotation on pulsation. Current models assume slow rota- 
tion with no effect on the pulsational geometry. The Coriolis 
and centrifugal accelerations become more important with 
increased rotation. Rotational effects can lead to differential 
rotation and, in extreme cases, the breakdown of spherical 



symmetry ( Townsend 2003 Lignieres et al. 2006 Reese et al 
2006). It is still not well understood to what extent rotation 



affects spectroscopic observations. Rotation can further hin- 
der frequency analysis as 7 Dor stars have rotational periods 
of the same order as their pulsational periods. 

This study focuses on one member, HD 135825, of the 
7 Dor group of stars to identify the frequencies and modes 
of the g-mode pulsations present and to look for any ob- 
servational evidence of hybrid pulsation and rotational ef- 



fects. Section [2] discusses the specific data acquisition and 
reduction undertaken. The analysis of the spectra follows in 
Sections |3] and |4] Specifically, Section |3] breaks down each 
of the frequency analysis processes. Section [4] then discusses 
the corresponding mode identifications to the frequencies de- 
tected. A discussion of results in context follows in Section 

El 



2 OBSERVATIONS AND DATA TREATMENT 

In total nearly 300 spectra of HD 135825 were collected from 
the 1-metre McLellan telescope at Mt John University Ob- 
servatory (MJUO) in Tekapo, New Zealand (170°27'.9 E, 
43°59'.2 S). The observatory is at an elevation of 1029 m. 
Spectra were obtained on the fibre-fed High Efficiency and 
Resolution Canterbury University Large Echelle Spectro- 
graph (HERCULES). HERC ULES operates over a w ave- 
length range of 3800-8000 A (iHearnshaw et al.' 2003 ) and 



the 4096x4096 pixel CCD installed in 2007 samples the en- 
tire free spectral range in one exposure. 

Multi-site data has also been taken for this star. A to- 
tal of 28 spectra were taken at on the Sadiford Cass Echelle 
Spectrometer on the 2.1 m telescope at McDonald Obser- 
vatory, USA. In addition 8 spectra were obtained on the 
SOPHIE spectrograph on the 193cm Telescope at the Ob- 
servatoire de Haute Provence, and 3 spectra on the McKel- 
lar Spectrograph on the 1.2-metre telescope at the Domin- 
ion Astrophysical Observatory. In order to incorporate these 
data a mean line profile with sufficient signal-to-noise for 
pulsation analysis must be constructed. The small size of 
the aforementioned datasets did not allow for a reliable mean 
profile with sufficient phase coverage to be produced. These 
data were therefore not incorporated in the present analysis. 

Data from MJUO were collected over a period of 18 
months from February 2009 to July 2010. These data were 
reduced using a MATLAB pipeline written by Dr. Duncan 
Wright. The pipeline performs the basic steps of flat field- 
ing from white-lamp observations, calculating a dispersion 
solution from thorium-lamp observations and outputting the 
data into a two-dimensional format. The observations were 
then processed in a pipeline written specifically for non- 
radial pulsation analysis. The processing stage corrects for 
barycentric motion, removes small differences between ob- 
servations, fits the continuum using a synthetic spectrum, 
then normalises and merges the orders for each spectroscopic 
observation. 

Each image was then cross-correlated using the delta- 
function method (Wright et al. (in preparation), building on 
Wright] ( |2008| ), [Wright et al.| ( |2007[ )) to create a representa- 
tive line-profile. The delta-function method relies on the fol- 
lowing assumptions of the spectra: spectral lines at different 
wavelengths have the same shape, lines of different species 
and different excitation potentials change in the same way 
as a result of the pulsation and all lines in a spectrum are 
varying in phase. The first two criteria were shown to hold 
to a good approximation in [Wright[ ( [2008[ ) and a check of 
lines of various equivalent width, species, excitation poten- 
tial, wavelength and showed them to be varying consistently 
in phase in a single spectrum for HD 135825. 

Time-series analysis and mode identification was done 
using the FAMI AS software ( Zima|2008 1 , applied as in Zima 



et al. (20061. This entailed an analysis of the moments 



I Balona 1986 Briquet & Aerts 2003 1 and pixel-by-pixel vari 



ations which identified the frequencies in the spectral lines. 
Specifically, FAMIAS uses the moment definition described 
in Aerts et al. ( 1992 1 and revised in Briquet & Aerts ( 2003 1. 



The software package SigSpec (Reegen 2007) was used 
as a validation of the frequency selection method. SigSpec 
performs a Fourier analysis of a two-dimensional dataset and 
selects frequencies based on their spectral significance. 

Mode identification was done using the Fourier Pa- 



rameter Fit method. Developed by Wolfgang Zima (Zima 



2009), this extension of the pixel- by-pixel frequency analy- 



sis method uses each identified frequency in the line profile 
variation and matches this to synthetic spectra with input 
stellar parameters for various modes until the best fit is ob- 
tained. 



3 FREQUENCY ANALYSIS 

HD 135825 is a relatively faint {V = 7.3), and hence com- 
paratively understudied, 7 D or candidate star. It has a mod- 
erate vaini of 38 ± 5kms~^ (De Cat et al. 20061. The star 



was confirmed multi-periodic in photometry by 



Eyer et al.| 



( 2002 1 and has an effective temperature of 7050 ± 90 K and 



log g of 4.39 ± 0.13, both from ([Bruntt et al.||2008[ ) and 
spectral type FO ( Kazarovets et al.|1999 |. A full abundance 
analysis has been done using high-resolution, broad wave- 
length range spectra by Bruntt et al. ( 2008 1 with a Metal- 



licity [M/H] = -1-0.13 ± 0.09 found. Full abundance re- 
sults are available in the on-line material. One pulsation pe- 
riod has previously been identified in photometry, 0.76053 d 
(1.31487 d"^), from the HIPPARCOS satellite data and one 
best fitting frequenc y from classification spectroscopy, 0.63 d 
(1.59 d"^), both by De Cat et al.| ( |2006] ). Note the main pe- 
riod from photometry was also suggested as the main period 
in De Cat et al. ( 2006 1 spectroscopy but was not recoverable 



due to the limited size of the radial velocity dataset. 

Each of the frequency variation tests (moments 0-3 and 
pixel-by-pixel) was used to construct periodograms for the 
star. Some methods have higher signal-to-noise and thus im- 
prove our ability to extract more frequencies. Amplitudes 
tabulated in this paper are scaled to ratios of /i to allow 
comparison between identification methods with differing 
absolute amplitudes. 

Representative line profiles are created by cross- 
correlating hundreds of spectral lines. These profiles of 
HD 135825 are shown in Figure IT] and are analysed using 
the pixel-by-pixel and zeroth to third moments. These tests 
cover all the independent frequencies. In addition to these 
tests, the software SigSpec was used to analyse all the two- 
dimensional data sets (all the moments) and also each indi- 
vidual pixel across the line profile, to confirm the frequencies 
found and their significance. 

The window function of the dataset was calculated to 
show any periodicities intrinsic to the data sampling. The 



window function for HD 135825 is plotted in Figure 2(a) As 



the dataset is from a single site the one-day aliasing is ob- 
vious and this can present a problem in the identification of 
frequencies. In this analysis, care has been taken to consider 
all one-day aliases and determine the true frequency. 




Velocity (kms ) 



Figure 1. Zero-point profile for the 291 observations of 
HD 135825 with the mean shown in red. The smooth mean shows 
that the variations in the profile are well sampled. The line profile 
data are provided as an on-line file that includes the Modified Ju- 
lian Date, velocity (on a relative scale) and the intensity at each 
of the 800 velocity sampling points for each profile. 



3.1 Pixel-by-Pixel 

The pixel-by-pixel method is an effective tool for examining 
the frequencies and their fits to the data in two-dimensions 
as this method performs an analysis of the zero-point line 
profile by considering the variation of each pixel. The cross- 
correlated line profile of HD 135825 contained 48 pixels. The 
pixel-by-pixel technique produced the lowest noise-level in 
the Fourier spectra. The Fourier spectra for the frequencies 



(fp) are given in Figure 2(b) and numerical results in Table 
[1] The phased variations of the zero-point profile are shown 
in Figure|3(a) and the fits for the four frequencies in Figures 
|3(b)|3(e)[ The profiles are smoothed and phased to show the 
variation over 1.4 cycles. Where a bin has more than one 
observation for a particular phase a signal-weighted mean 
has been used. These plots show the increased matching to 
the data shape for each successive additional Frequency with 
the first frequency dominating the shape and the fit. 

Table [l] column five shows the amount of variation ex- 
plained by each frequency combination calculated from the 
residuals of the fits to the data at each pre-whitening stage. 
The percentage of the variation explained in the pixel-by- 
pixel method is generally lower due to several factors. Firstly 
we are considering the motion of all 48 pixels in the line 
profile, only about 15 of which are strongly varying (see the 



amplitude variation profiles as in Figure 9(a) I. This means 



we are considering high and low signal-to noise regions of 
the line profile. Another effect of considering all pixels in 
the line profile is the increased precision required to match 
all the individual pixels. Frequencies, amplitudes and phases 
must be much more precise for an accurate sum of pixel fits 
to account for the movement in the line profile. Despite this 
the Fourier spectra and the reduction of residuals in each 
pre-whitening stage leaves us confident in the identification 
of these frequencies. With the fit of four frequencies, some 
29% of the variability in the line profile remains. Much of 
the remaining variation is likely due to further un-identified 
frequencies at much lower amplitudes. The signal-to-noise 



limit of our data (see the red line in Figure 2(b) I restricts 
the extraction of further frequencies. 



Table 1. Frequencies from the pixel-by-pixel analysis of 
HD 135825 from MJUO. 



ID 


Frequency 


Amplitude 


Phase 


Cumulative 




id-') 


(scaled to /pi) 




Variation Explained 


/pi 


1.3149 


1 


0.1033 


49% 


,fp2 


0.2901 


0.4786 


0.4199 


58% 


/p3 


1.4045 


0.4717 


0.0368 


66% 


/p4 


1.8830 


0.3874 


0.4058 


71% 




Frequency (d ^) 
(a) 




3 4 5 

Frequency (d" 



(b) 



Figure 2. The Fourier spectra showing the amplitude of fre- 
quencies detected in the d ata. | (a) [ Window function for the 291 
observations of HD 135825. |(b)| frcquencies found using the pixel- 
by-pixel technique and successive pre-whitening. The smooth line 
shows the significance level taken to be the limit of detection. 
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(a) The dataset from MJUO phased to /pi = 1.3149 d^^ 
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(b) Fit with /pi (1.3149d-i) (c) Fit with /p2 (0.2901d'i) 
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(d) Fit with /p3 (1.4045d-i) (e) Fit with /p4 (1.8830d-i) 

Figure 3. Residual (line profile minus the mean line profile) 
phased line profiles and fits to the profiles usin g fr equencies (in 
d^^) identified with the pixel-by-pixel method. |(a)| Real dataset 
and |(b)| to |(e)| are the fits to the profiles with the addition of suc- 
cessive frequencies and modes constructed as described in Section 



3.2 Zeroth Moment (Equivalent Width) 

Variations in tlie equivalent widtii of the line profile are 
not expected unless there are significant temperature vari- 
ations, which makes the analysis much more complex. 
For HD 135825 only one frequency could be extracted 
(1.3147 d~^), with a low-amplitude Fourier peak above the 
noise base. The lack of frequencies in the zeroth moment 
indicates that there are only small temperature variations 
centred around the primary frequency and that this star is 
a good candidate for spectroscopic pulsation analysis. 



3.3 First Moment (Radial Velocity) 

The radial velocity measure is a good frequency analysis 
tool as it is very sensitive to any low-degree modes. First mo- 
ment results are in good agreement with the results from the 
pixel-by-pixel search. The Fourier spectra for the frequencies 
(/a) are given in Figures 4(a)|4(c) and numerical results are 



shown in Table [2] Due to an uncertainty in the selection of 
the Fourier peaks, three paths of frequency were selected, 
A, B and C. The first frequency identified, 0.3122 d~^, was 
a 1-day alias of the main frequency discovered in all other 
methods. The second peak, 1.3150 d~^, had an amplitude 
almost identical to the first peak and the latter was taken as 
the true first frequency for analysis. Increasing amplitudes of 
Fourier spectra is a known effect of higher noise levels in the 
to 1 d~^ region of the Fourier spectrum, which can push 
an alias peak of a frequency higher than the real frequency. 
There is an ambiguity in the selection of the second and 
third frequencies. There are two possibilities for the second 
frequency, either 0.2412 d"^ or 0.2902 d"\ The 0.2902 d"^ 
frequency was chosen as the other is not consistent with the 
pixel- by-pixel and SigSpec results. The 0.2902 d~^ peak is 
the second highest peak in the pixel-by-pixel method. Pre- 
whitening using the 0.2902 d~^ frequency removed this peak 
entirely from the Fourier spectrum in both methods so it 
was regarded as the best choice for fa2- The 0.2412 d^^ 
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Figure 4. The Fourier spectra showing frequencies found using 
the 1** moment technique with successive pre-whitening. The red 
hne shows the significance level taken to be the limit of detection. 
Path C was taken to be the best selection of frequencies due to the 
improvement in the percentage variation and overall frequency 
amplitude. 



frequency was still considered in the analysis of Path A. 
The third frequency is more complex as it first appears as 
1.4045 d~^ in the pixel- by-pixel method but then a very 
close frequency, 1.4176 d~^, appears using the radial velocity 
method. It is worth noting that in each of the Fourier spectra 
the other frequency appears strongly as an additional peak 
and that choosing the same frequencies as found using the 
pixel- by-pixel method gives the same frequency set. When 
1.8831 d~^ was found using this method it had a higher 
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Figure 5. Radial velocity measurements and fits in a section of 
the HD 135825 data. The crosses show the observed radial velocity 
and the solid and dashed lines show the fits using Path A and Path 
C frequencies respectively. 



amplitude relative to fai than /a4 given by choosing just 
the highest peaks. This confirms our choice of Path C. 

The radial velocity data is plotted over a selected range 
of Julian date in Figure [S] This Figure gives an indication 
of the difference in fit between Path A and Path C. 

Even after the removal of four frequencies there still ap- 
pear to be peaks and aliases in the Fourier spectrum that 
have amplitudes that are too low to significantly distinguish 
them from the noise. This could indicate further frequen- 
cies are present or that one or more frequencies have been 
slightly misidentified. 

Table|2]shows the amount of variation explained by each 
frequency combination. 



3.4 Second Moment (Variance) 

The second moment refiects variations in the width of the 
line profile and is particularly sensitive to even values of 
m. In general we expect the first and second moments to 
be complementary and together build a frequency list sim- 
ilar to the pixel-by-pixel findings. In the analysis the first 
frequency identified was the usual 1.3149 d~^ and the sec- 
ond frequency that was recovered was 0.4398 d~^. Another 
high peak was seen at the previously identified frequency 
0.4043 d^'^ so both options were studied as a possible fre- 
quency path. Frequencies (fb) are tabulated for both path 
D and E in Table |2] Path E found a higher relative ampli- 
tude in fb3 than Path D. In the fb4 spectrum a 1-day alias 
of the previously identified 0.290 d~^ is seen as a high peak 
in both Fourier spectra. Path E most closely resembles fre- 
quencies identified in the pixel-by-pixel method. Tableplalso 
shows the amount of variation explained by each frequency 
combination. 



3.5 Third Moment (Skewness) 

The third moment frequency analysis is known to produce 
similar frequency sets to the first moment method. It de- 



Table 2. Frequencies and amplitudes, relative to fia, fib, fic identified using the variation of the I'**, 2 and 3^ moments respectively. 
Note (n) identifies the rank of the frequency peak (high-low) in the case where the highest peak was not chosen. Paths indicated in bold 
indicate the best choice of frequencies. See text for further explanations of the frequency paths. 









First Moment 






ID 




Frequency (d ^) 




Amplitude 


Variation Explained 




Path A 


Path B Path C 


Path A 


Path B Path C 


Path A 


Path B Path C 


fal 
fa2 
faS 
fa4 


0.2412 
1.4124 
0.9788 


1.3150 

0.2902(2) 
1.4149 0.4040(7) 
2.4105 1.8831 


0.5351 
0.3931 
0.3206 


1 

0.5278 
0.4092 0.3388 
0.3085 0.3761 


69% 

76% 
83% 


53% 

67% 

77% 75% 
84% 83% 



Second Moment 





Path D Path E 


Path D 


Path E 


Path D Path E 


.fM 


1.3149 




1 


63% 


.fh2 


0.4398 0.4043(2) 


0.3675 


0.3421 


70% 70% 


fb3 


1.8829 1.8829 


0.3401 


0.3519 


78% 78% 


hi 


0.2430 0.7125 


0.2454 


0.2666 


82% 81% 









Third Moment 








PathF 


Path G Path H 


PathF 


Path G Path H 


Path F 


Path G Path H 


.fcl 




1.3150 




1 




58% 


.fc2 


0.4126 


0.2903(3) 


0.4330 


0.3948 


69% 


67% 


.fc3 


0.2551 


0.4150 0.4040(8) 


0.3417 


0.3437 0.3142 


76% 


77% 75% 


/c4 


3.4104 


2.4104 1.8831 


0.2546 


0.2807 0.3429 


83% 


84% 83% 



scribes the variation in the skewness of the line profile. In 
this study the third moment also produced the most com- 
plicated frequency spectra with three paths (Path F,G and 
H) or frequency (/c) sets being identified (see TablelSJ. This 
complication was due to the ambiguity in peaks in the iden- 
tification of fc2 and /c3. In this study Path H is chosen as 
that which identifies the same frequencies as the pixel-by- 
pixel method and had the highest relative amplitude of /c4. 
The frequency 0.41 d~^ was found preferentially in /c2 and 
/c3 but vanishes if 0.405 d~^ is chosen. The 0.405 d~^ is still 
present in Path F and G as it (or a 1-day alias) occurs as 
/c4 and /c3 respectively. This supports the choice of Path H 
as the best set of frequencies. Table [2] shows the amount of 
variation explained by each extracted frequency set. 




3.6 SigSpec Results 



The software package SigSpec ( |Reegen||2007[ ) was used to 
analyse the two-dimensional moments for frequency identifi- 
cation and significance. SigSpec searches the discrete Fourier 
transform for frequencies and assigns spectral significance. 
Spectral significance is a measure of false-alarm probability 
for a given frequency. This probability refers to the likeli- 
hood that random noise in the time domain could generate 
a peak in the Fourier amplitude of similar size as that of 
the data itself. For example if the risk a noise peak ap- 
pears at this amplitude is 1:100 000 then the false-alarm 
probability is 0.00001. This value is equivalent to a spec- 
tral significance of 5.0 ( Reegen|2007 l. The default threshold 
value for determining frequencies in SigSpec is set at 5.46 
(equivalent signal-to-noise = 4) but for the purposes of this 
analysis we took a much larger threshold of 10. This caution 
reflects the requirement of precise frequency determination 
for accurate mode identification. For the three-dimensional 
pixel- by-pixel dataset, all 48 pixels were individually anal- 
ysed for frequency identification using SigSpec. The results 
for each moment are given in Table |3] for all frequencies 



Figure 6. The 2D peridogram (left) showing frequencies found 
across each pixel (lighter colours for stronger peaks). The peaks 
at 1.3 d~^and 1-day aliases of this frequency are clearly visible. 
The mean frequencies of the 48 pixels is plotted on the right. 



(/s) with spectral significance greater than 10. Many of the 
frequencies found with significances between 10-15 are not 
reproduced consistently in other methods of frequency de- 
termination and these may not represent true frequencies in 
the stellar pulsation. 

Figure [6] shows the Fourier frequency identification 
across all 48 pixels and shows the three-bump structure 
of the line-profile variation and the dominant frequency at 
1.315 d"'^ and 1-day aliases. 



3.7 Frequency Aliases and Combinations 

As this analysis produced a significant set of spectroscopic 
frequencies for a 7 Dor star, a full study of the Fourier spec- 



Table 3. Frequencies found using SigSpec above the significance 
threshold of 10. Frequencies that appear in multiple analysis 
methods are indicated in bold. 





1"* mom 


Sig 


2"'' mom 


Sig 


3'''' mom 


Sig 


fsl 


1.3151 


33.8 


1.3151 


40.4 


1.3151 


37.9 


.fs2 


0.2906 


19.1 


1.4425 


13.2 


0.4125 


16.1 


/s3 


0.4149 


17.0 


1.8834 


15.3 


0.2541 


14.2 


.fs4 


1.9262 


15.4 


0.2435 


11.8 


2.4076 


12.5 


fsS 


1.4457 


15.2 






1.9685 


12.7 


/s6 


2.5257 


11.9 






2.4720 


10.5 


fs7 


2.1490 


11.6 






2.1855 


12.3 


fs8 










2.6724 


10.5 


fs9 










0.0283 


11.7 


falO 










0.2013 


11.3 



tra was undertaken. The presence of some combinations for 
particular frequencies led us to investigate aliasing patterns 
and possible combinations of true frequencies which may 
show up in our data. This was done using the frequencies 
found using the pixel-by-pixel method although other inde- 
pendent frequencies found in the moments were considered. 

The 1 cycle-per-day alias pattern is obvious in our data 
(Figure |2(a)[ ) due to observations being taken from a sin- 
gle site. Aliasing however does not seem to add significant 
uncertainties to our frequency identification as the aliases 
have significantly lower power than the identified frequency. 
The case where this is most evident is the identification of 
/p3. Mode identification tests of both frequencies showed 
1.405 d~^ was a better fit to the line profile variation than 
0.405 d~^. This was not the case for fp2 where the 0.290 d~^ 
frequency was found to be a better match to the data, de- 
spite being in the longer range of periods for a 7 Dor star. 

Combinations of frequencies are known to occur in- 
herently in the star. The phenomenon has previously been 
taken into account in the analysis of photometry of 5 Scuti 
stars such as in Breger et al. (20051. For our set of four 



identified frequencies a grid search of combinations includ- 
ing addition and subtraction of all frequencies, their multi- 
ples and their one-day aliases was undertaken. It was found 
that there is possibly a link between fp2, fps and fp4 as 
(/p2 + 1) — fp3 ~ fp4 — 1- It was also apparent that the first 
identified frequency for /s in both the first and third moment 
is not found indicating it could be an independent frequency. 
The frequency /2 found first in the second moment is half 
that of identified frequency fp4 (excluding one-day aliasing). 




Frequency d ^ 



Figure 7. The Fourier spectrum of the synthetic data with only 
one implanted frequency of /pi = 1.3150d~^ at the observational 
spacing. 
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Figure 8. The four Fourier spectra for the identified frequencies 
in the synthetic reproduction (left) and in the real pixel-by-pixel 
dataset on the right. 



they are artefact frequencies as only the above four frequen- 
cies inserted should be recovered in the Fourier spectrum. 



3.9 Frequency Results 

A summary of the frequencies found by all methods is pre- 
sented in Table HI An error estimate for each frequency 
((t(/)) is given using the definition of Kallinger et al. (20081 
who propose: 



3.8 Synthetic Frequency Identification 

As an additional test of the validity of the derived frequen- 
cies a synthetic data set was constructed using FAMIAS 
( |Zima|pb08[ |. The data used the same time spacings as the 
real observations and then the four frequencies, their ampli- 
tudes, phases and modes (as determined in Section |4| were 
entered along with the known parameters of the star and 
representative spectral line. The results for the derivation of 
the first frequency alone in the data and the first four fre- 
quencies together are given in Figures [7] and [8] respectively. 
All frequencies implanted in the line profiles were able to be 
extracted, although /p4 appeared in the Fourier spectra as a 
higher peak than /p3. The frequencies 2.80 d~^ and 1.11 d^^ 
are found as the fifth and sixth frequencies which indicates 



<^(/) = 



T* \/sig{a) 



(1) 



where T is the time base of observations in days and sig{a) is 
the spectral significance of the frequency from SigSpec. This 
equation provides a low uncertainty on our frequencies due 
to the long time base of observations, but should be regarded 
as a guide only as the authors feel this uncertainty under- 
represents the errors present in the frequency identification. 
A further least-squares fit of the frequencies to the one- 
dimensional data was performed to get a second estimate 
on frequency errors by examination of the 95% confidence 
bounds. The bounds gave errors on the order of ±0.00002 
d~^. These errors were disregarded and the previously found 
SigSpec uncertainties adopted from this point as they are 
an order of magnitude larger than the confidence fits and it 



is prudent to remain cautious about the precision of these 
results. 

Four frequencies are confirmed in this study: 1.3150 ± 
0.0003 d-\ 0.2902 ± 0.0004 d'S 1.4045 ± 0.0005 d"^ and 
1.8829 ± 0.0005 d~^ (Table |4| These are the frequencies 
found originally in the pixel-by-pixel measurements and con- 



firmed in the 1" 



moment analysis with uncertain- 



ties calculated from significances in the 1" moment. Ad- 
ditional frequencies that may be present in the data are 
1.44 d"\ 0.25 d"^ and 1.415 d"^ These frequencies had ei- 
ther poor signal-to-noise or misshapen amplitude variation 
profiles which limited any further mode-identification. 

From the identification of multiple frequencies in the 
0.3-3 d~^ range we can confirm HD 135825 is a true 7 Dor 
star. After the removal of four frequencies in all the above 
methods there remained some indication of a periodic signal. 
It is likely there are further 7 Dor frequencies present in the 
data below our detection limit with this dataset. There do 
not appear to be any higher frequencies suggestive of a 5 
Scuti hybrid. 



4 MODE IDENTIFICATION 

The first calculation performed was a fit to the zero-point 
profile using the parameters vsini, equivalent width and ve- 
locity offset. The latter is used to bring the observational 
data onto a consistent scale. These values were then used 
to refine future searches in parameter space. To continue 
with mode identification of the pulsation frequencies, the 
three Fourier parameters (zero-point, amplitude and phase 
variations across the cross-correlated line profiles) were cal- 
culated. Using the Fourier Parameter Fit Method (jZima 



2009), the modes of each of the final frequencies in Table 4 



were found independently, as well as the combination of the 
four frequencies. This was done using the frequency analysis 
package in FAMIAS (|Zima||2008|). The fits for each mode 



were performed on the amplitude and phase distributions 
and a x^ calculated. The pulsation parameter space that was 
searched (using a genetic algorithm optimisation to search 
parameter space) is outlined in Table [5] 

Figures 9(a)|9(dy show the Fourier parameters and their 
fits to the amplitude and phase across the line profile for the 
four frequencies. Note that FAMIAS sometimes generates 
phases (see Figures 9(a)|9(b)|9(d) \ that are ±1 cycle (±27r) 
displaced with respect to the observations. The mode fits are 
best for the even frequencies, /2 (=0.2902 d"^; (l,m) = (2,-2)) 
and /4 (=1.8829 d'^; (1,1)) that have x^ = 3.1 and 2.5 re- 
spectively, whereas /i (=1.3150 d'^ (1,1)) and /g (=1.4045 
d"^ (4,0)) have x^ of 12.9 and 10.6 respectively. The most 
significant deviations are in the /a fits where the amplitude 
and phase are not fitted well, albeit with lower pulsational 
amplitudes than /i and /2. 

The final parameters of the fit to each of the modes in- 
dividually, the x^ value and stellar parameters are given in 
Table \g\ The mode with the lowest x^ fit was selected as 
the best mode. This was an unambiguous selection as there 
was only one mode for each frequency with a significantly 
lower x^- The x^ o f each mode fit to the final frequencies 
is given in Figures 10(a)|lO(d) The individual x^ values 



Table 5. Stellar parameters used in the mode identification. 
Fixed values for Teff, [M/H] and log g were taken from |Bruntt| 
|et al.| | [2008| and vsini values from pe Cat et al.| ( |2006^ . Values for 
the other parameters were taken as reasonable limits for 7 Dor 
stars. 



Property 



Fixed Value Min Max Step 



Radius(solar units) 




1 


5 


0.1 


Mass (solar units) 




0.5 


5 


0.01 


Temperature (K) 


7050 








Metallicity [M/H] 


0.13 








log 9 


4.39 








Inclination (°) 







90 


1 


tisinj (kms^-"-) 




30 


50 


1 



Table 6. Results of mode identification for /i alone, all four 
frequencies individually after a least-squares fit is applied (Isf), 
and all four frequencies simultaneously (sini). 



frcq 


Mode ID 


x' 


Inclination 

(°) 


Vel. amp 
(kms-i) 


Phase 


/i 


(1,1) 


7.6 


87 


0.78 


0.64 


/i Isf 


(1,1) 


12.9 


64 


0.78 


0.64 


/2lsf 


(2,-2) 


3.1 


55 


0.65 


0.43 


/alsf 


(4,0) 


10.6 


88 


0.65 


0.40 


/4lsf 


(1,1) 


2.5 


51 


0.65 


0.88 


/i sim 


(1,1) 


12.7 


87 


0.80 


0.63 


/2 sim 


(2,-2) 






0.50 


0.41 


/s sim 


(4,0) 






0.67 


0.40 


fi sim 


(1,1) 






0.69 


0.91 



fit constraints. Couple this with the poorly modelled asym- 
metry in the line profile and the x^ can be higher than is 
expected with a good fit such as for /i. When all modes 
were combined for a simultaneous fit a reduced x^ of 12.7 
was achieved with the stellar variables and line parameters 
in TablelS] The other best possible fits to the data are shown 
in TableTTl Although some alternate fits to /a have a sim- 
ilar x'^, the shape of these other fits do not show the four 
bumps as in the line profile variation. The second and third 
best fit modes have three and two bumps respectively in the 
amplitude variation profiles and do not display four distinct 
changes in the phase and thus we are confident the (4,0) 
identification is the best match. 

The best individual values for the parameters are given 
by the best fits which are to the second and fourth frequen- 
cies. The values for the first frequency also merit consid- 
eration given its dominant contribution to the line-profile 
variation. The range of values found for the inclination of 
the star show it is an uncertain result of the search, but 
values between 50° — 90° are reasonable and an inclination 
close to 90° may explain the number of detected frequen- 
cies, as we observe the star edge-on. This means we see the 
full surface of the star as it rotates. The value for v sini was 
measured from the line profile to be 38 ± 5 kms"'^ (De Cat 
et al.|2006| ). 



can be used as a guide to the goodness of fit but frequen- 
cies with higher amplitudes and snrall errors have tighter 



4.1 Rotation and Pulsation Parameters 

FAMIAS has a tool to provide a check for the validity of the 
mode identification method based on the rotational param- 
eters of the star. The stellar rotational parameters of equa- 



Table 4. Frequencies (in d~^) identified using pixel- by-pixel (pbp) and 0*'' to 3'''^ moment methods. Note (a) denotes where a one- 
cycle-per-day alias was found as the highest peak and (#) identifies frequencies which did not have the highest peak but were sequenced 
according to the number in brackets and chosen as a better result. 





pbp 


0th 


1st 


2nd 


3rd 


final (d-i) 


.fl 


1.3f50 


1.3147 


1.3122 (a) 


1.3149 


1.3150 


1.3150 ± 0.0003 


.f2 


0.2902 




0.2902 (2) 


0.4043 (2) 


0.2903 (3) 


0.2902 ± 0.0004 


h 


1.4045 




0.4039 


1.8829 


0.4040 (8) 


1.4045 ± 0.0005 


U 


1.8829 




1.8831 




1.8831 


1.8829 ± 0.0005 
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(a) /i = 1.3150 d-i fit with (1,1) mode, x^ 



60 
12.9. 





(b) /2 = 0.2902 d-i fit with (2,-2) mode, x^ = 3.1. 
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Figure 9. Mode identification results. Amplitude variation and phase across the line profile (solid) of the four identified frequencies with 
the lowest x^ fits (dashed). A maximum error bar for the fits is given in each amplitude variation plot. Errors on the data and phases 
are too small to display. 



Table 7. Best fits for modes identified for HD 135825 



/l 



h 



h 



(1,1) (2,-2) 



(4,0) (1,1) 10.6 


(1,1) 


12.7 


(2,-2) 


13.9 


(5,4) 


14.9 


(3,-1) 


15.5 


(3,-3) 


15.6 


(4,3) 


15.8 



(vcrit) , critical vsini (vsadcrit) and critical minimum incli- 
nation (i) are calculated from inputs of mass, radius, wsini 
and inclination using estimates for mass and radius typical 
of 7 Dor stars (1.5 M© and 1.6 R©) and usini and inclina- 
tion from the mode identification (39.7 kms~^ and 87°). The 
results are tabulated below: 



torial rotational velocity (vrot; ^rot = 27rR/rot), rotational 
period (Trot), rotational frequency (frot), critical velocity 



Table 8. Stellar rotational frequencies. Final frequencies, their 
co-rotational values (fco-rot), horizontal-to-vertical amplitude 
ratio (k, see equation[2]| and ratio of rotational frequency (frot) 
to each co-rotational frequency. 



Frequency 


Frequency 


J CO — rot 


K 


Jrot 1 Jco — rot 


ID 


d-i 


d-i 






h 


1.3150 


0.8241 


40 


0.6 


h 


0.2902 


1.2720 


17 


0.4 


h 


1.4045 


1.4045 


14 


0.3 


h 


1.8829 


1.3920 


14 


0.4 



^rot 
-L rot 

Jrot 
^crit 

"^crit 



39.75 kms"^ 
2.03 d 
0.49 d-i 
423 kms"^ 
423 kms"^ 
5.4° 



From the above table we can see that HD 135825 is not near 
the break-up limit of the star and is naturally rotating with 
a period of around two days. The ratio of horizontal to ver- 
tical amplitude of the pulsation (k) , is calculated using the 
equation 



GM 



(2) 



where G is the gravitational constant, R and M are the stel- 
lar mass and radius and u is the angular frequency. For the 
various frequencies k, was much greater than 1.0 which in- 
dicates low-frequency g-modes as we expect for 7 Dor stars. 
We can also use the determined inclination to transform 
our observed frequencies into the co-rotating frequencies 
ifco-rot) required to model the star {fco-rot = fobs — mfrot)- 
These frequencies, along with the horizontal-to-vertical am- 
plitude ratios calculated using Equation[2J are given in Table 
Is] All of the frequencies lie now in the range accepted for 
7 Dor stars which strengthens our confidence in having ob- 
tained correct frequency identifications, particularly for /2 
which was unusually low. We also show the ratio of the rota- 
tional frequency to the co-rotating frequency which can give 
an indication to the validity of the mode identification mod- 
els being applied by FAMIAS. In general ratios less than 
0.5 are considered to lie well within the parameter space 
which FAMIAS can produce reliable results, however this is 



dependent on the value of m considered. Townsend ( 2003 1 



demonstrated that prograde modes (m > 0) are distorted 
less from increasing Coriolis force, increasing the limits of 



reliability. See Wright et al. (20111 for further discussion. 



Note the results obtained in this section must be re- 
garded as approximate as we have estimated the stellar mass 
and radius. 



5 DISCUSSION 

The dataset in this study provides the largest spectroscopic, 
and indeed largest observational, dataset for HD 135825 to 
date. We are now in a position to examine our results in the 
context of previous work on this star. Multi-periodicy was 



on page 209 HD 135825 (HIP 74825) is incorrectly written 
as HD 135828) but no frequencies were published. One pul- 
sation period was identi fied in HIPPARCOS photometry, 
0.76053 d (1.31487 d"^) ( |De Cat et al]|2006[ ). In the same 
paper 17 spectra were analysed and one period peak identi- 
fied at 0.63 d (1.59 d~^), but this was treated with caution 
due to the limited spectra available at the time. No mode- 
identification has previously been published. As mentioned 
above, the wsini value of 39.7 kms ~^ matches we l l with pre- 
viously published of 38 ± 5 kms" 



De Cat et al. 



(20061 



The first frequency extracted in this study matches that 
found in HIPPARCOS data, but the 1.59 d^^ frequency was 
not found. 

The high numbers of single-site data can be used to 
extract multiple frequencies from a g-mode pulsator. It is 
worth noting that much of the ambiguity in the path selec- 
tion and hence frequency identification in the i^'-3'"'* mo- 
ments could be removed using multi-site data. It is also nec- 
essary that any multi-site data is acquired such that it covers 
the phase range of the frequencies independently. Figure [5] 
shows how multi-site data could distinguish clearly the bet- 
ter frequency fits by providing data points where the fits 
differ. Despite the increased difficulty in extracting the pul- 
sations from single-site data, the authors are confident that 
the frequency and mode identifications given in this research 
are robust, given that similar results are obtained from in- 
dependent analysis techniques. 

We can now begin to remark on the size and quality of 
spectra required to undertake geometric pulsation analysis 
in spectroscopy. The dataset produced high signal-to-noise 
frequencies and clear mode identifications using 291 high- 
resolution spectra taken over 18 months. It is clear that 
this is a demanding standard to study all 7 Dor stars but 
the significant increase in precision and ability to success- 
fully identify multiple frequencies warrants this approach. 
It is noted that not all stars will produce similar results 



with large datasets. For example, HD 40745 ( Maisonneuve 



et al.|[2011j has complicated frequencies and modes even us- 



ing more than 400 spectra of a comparable quality to this 
study. Even more spectra (nearly 700 from multiple sites) 



were obtained by Uytterhoeven et al. ( 2008 1 on the star 



HD 49434 and up to six frequencies were found in each of 
the spectroscopic methods although only some frequencies 
were found in multiple methods. 

Single-line analysis of the frequencies in HD 135825 pro- 
duced the same /i in phase for all lines tested. The assump- 
tion that all spectral lines move in phase has been found 



true for all 7 Dor stars published to date (see Maisonneuve 



etal. ( 2011 1 , Wright ( 2008 1 for more examples) . This demon- 



first confirmed in photometry by Eyer et al. ( 2002 \ (Note 



strates that these stars are suitable for application of the 
cross-correlation method and this technique is recommended 
to obtain sufficiently high signal-to-noise to determine mul- 
tiple frequencies in spectroscopic data. 

As observational asteroseismology of 7 Dor stars moves 
beyond the classification phase into producing more indi- 
vidual identifications, it seems likely most stars will require 
similarly large, or larger datasets as this study. Stars with 
temperature variations and other periodic structure, such 
as tidal effects from binaries, will provide further challenges 
and also require massive datasets. 

It is clear that high-precision spectroscopic mode iden- 
tification is dependent on the availability of high-precision 
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(a) x^ values ior I = - 3 for /i = 1.3150d-i. 
Best fit value is for {l,m) = (1,1). 



(b) x^ values for « = - 3 for /2 = 0.2902d-i. 
Best fit value is for {l,m) = (2,-2). 
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(c) x^ values for i = - 4 for /i = 1.4045d-l. 
Best fit value is for {l,m) = (4,0). 



(d) x^ values for « = - 3 for /i = 1.8829d-i. 
Best fit value is for (l,m) = (1,1). 



Figure 10. Lowest x^ values for each possible (Lm) combination for the final identified frequencies. The best fit x^ is identified in bold. 



spectral abundance analysis and modelling. The production 
of the best mode identification is reliant on the availability 
of Te// and log g measurements and also benefits from good 
estimates of stellar masses, radii and inclination. The pro- 
duction of co-rotating frequencies is crucially dependent on 
precise values of inclination and radius, which are not well 
constrained in the current mode-identification method. 

Even small numbers of full mode-identified frequencies 
can be used to place constraints on stellar models, particu- 
larly on the conditions required to produce such a selection 
of excited modes. Future work is to take the mode identifica- 
tions of HD 135825 from this paper and use them in complex 
theoretical models (e.g. |Townsend|2003[ ) to further constrain 
stellar parameters and to study the extent of the core. This 
could give us important insights into the evolutionary past 
and future of the star. Feedback from such models will ad- 
ditionally inform our mode identification methods and lead 
us ever closer to understanding these information-rich 7 Dor 
stars. 
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